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This paper examines several aspects of the discrete matrix model of population transition. Certain 
appropriate applications of matrix theory and exploitation of the specific form of the model should 
serve to enhance its already well-developed status. The aspects dealt with include (1) a simplification 
of the Perron-Frobenius theory; (2) row and column sum bounds on maximal eigenvalues; (3) relations 
between oscillations in a population and the remaining eigenvalues; (4) implications of stability for 
the transition matrix; and (5) relations between characteristic quantities of a stable population. 
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1. Introduction 

The purpose of this paper is to examine certain aspects of the matrix model of population 
dynamics. We begin with a brief discussion of the model. More details may be found in Keyfitz 
(1968). 

In the discrete model of population, the population at time t is represented by a vector, v{t). 
The fcth component of this vector, vi(t), is the number of people at time t between the ages of 
(i — I) At and iAt where At is a fixed time period, usually one or five years. 

The following matrix A= (atj) acts onv(t) togivei;(£+ 1), see Keyfitz (1968): 



a u a l2 
a n 
a 32 



a Ul 







a „ 







(1) 



where a xu . . . , a Xn are age-specific fertility rates and a 2 i, . . . , a m ,m-i are survivorship ratios, 
both of which are based on the time period At. In this paper we assume these rates remain constant. 
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A matrix, X, with nonnegative entries is called reducible if there exist a permutation matrix P 
such that 



P l XP 



Y 
W 



(2) 



when Y and Z are square matrices. Matrices which are not reducible are called irreducible. The 
matrix (1) is reducible. However if we restrict ourselves to the first n rows and columns we get the 
following irreducible matrix 



A = 



an . . 


• 


a ln 


«21 







a 32 













a„,„-t 



(3) 



where am is the last nonzero fertility rate. We also restrict v(t) to the first n components. This 
restricted model still contains all the essential information, see Parlett (1970). 
The most important theorem for this population model is 

THEOREM 1: (Perron-Frobenius): Let A be an irreducible nonnegative (entry-wise) matrix. Then A 
has a simple positive eigenvalue r called the Perron root with corresponding eigenspace of dimen- 
sion 1 and such that the absolute value of all other eigenvalues is less than or equal to r. If A has 
exactly h roots with absolute value r then these roots are: r, re l27r/h , re l47r/h , . . ., re l(h_i)27r/h . 
Furthermore, an eigenvector associated with r can be chosen to have positive components. 

The proof of this theorem uses the Brouwer fixed point theorem. However, parts of this proof 
can be simplified for the matrix^ in (3). The characteristic polynomial of A is 



f(\) = k" — auk"- 1 — a X2 a 2X \ n - 2 — a l3 a 32 a 21 \ n - 3 . . . — am#n,n-ian-i,n-2 . . • a 2 \. 



(4) 



Since /(0) < and/(*) > for sufficiently large*, there is a positive root r, that is/(r) =0. 
To prove that/(*) has no root with greater modulus than r we need the following theorem due to 
Rouche see Nehari (1961), p. 138. 
THEOREM 2 : If gi (z) and g 2 (z) are polynomials over the complex numbers such that 

| gl (z)|>|g 2 (z)| 

everywhere on a circle C, then the polynomial gi(z) -h g 2 (z) has the same number of roots, counting 
multiplicity, inside C as gi(z). 

It is clear from (4) that/ (R ) > for all R > r. Hence 



R n > a X iR n ~ A + a 21 ai 2 R n ~ 2 + a 21 as 2 a 13 R n ~ 3 + . . . + a 21 a 32 . . . a n ,n-\a\,n 



(5) 



for all R > r. Let gi{z) = z n and g 2 (z) =f(z) — z n . Let R > r, then if z is on the circle of radius R 
about the origin, we get: 
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1^2 (z) | = |— (i\\z" ' — a 2 ia V2 z"- 2 — . . . — a 2 i«32 • . - a n ,n-\(i\n\ 

«s flnfi"-' + a 2 ia 12 R n ~ 2 + . . . H-a 2 ia 3 2. . . a R , n -ia ln (6) 

< /?" by (5). 

It follows that all the roots of/(z) lie within the circle |z| = R. Since R can be made arbitrarily 
close to r, no eigenvalue can exceed r in absolute value. We also point out that r lies between the 
minimum and maximum row sums and also between the minimum and maximum column sums 
of (3). Since the row and column sums of (3) are of a very simple nature, in specific cases this 
may give fairly tight bounds on r. For example, in a human population the survivorship ratios a,,/-i, 
i = 2, . . . , n may fall within a small range, and if total fertility, the sum ci\\-\- . . . + a ]M also falls 
in that range, then r must fall in that range. For the proof of these and references to further inequal- 
ities see Gantmacher (1959). 

In general even better bounds on r can be attained from the a\j at the price of greater computa- 
tion, see Marcus and Mine (1964). Perhaps the best of these applicable to an irreducible matrix 
of the form (3) is the following. 

THEOREM 3: Let s be the minimum row sum and S the maximum row sum of A. Let a be the sum 
of all entires of A divided by n, and let k be the smallest positive off diagonal entry of A. Set 
e=(A/S)"- 1 . Then 

s + e(cr-s) ^ r^ S-e {S - a) . 

We make one more definition which will be used later. An irreducible nonnegative matrix is 
called primitive if its Perron root is strictly larger than the absolute value of all other character- 
istic roots of A. A matrix which is not primitive is called imprimitive. A detailed discussion of 
primitivity is contained in Sykes (1969). 

2. Oscillations in Population 

It is well-known that under constant fertility and mortality the age specific population vector 
converges to a multiple of the eigenvector corresponding to the Perron root, provided that the 
population matrix is primitive. In this section we investigate what is happening to the population 
vector as it approaches its stable value. We will show that the population vector oscillates much 
like a damped pendulum, eventually coming to "rest" at the stable population vector. 

Bernadelli (1941) discusses the importance and causes of these waves. 

Let \i = r, A.2> • • • , k n be the eigenvalues of the matrix A given in (3), numbered so that 
\ x 55 |X 2 | ^ . . . ^ |X n |. Let r, = |X,| and let 0,- be the angle \% forms with the x axis, — 77 < 0; ^ tt. 
Assume for now that the eigenvalues of A are distinct and that if r, = rj then i =j or r, is the complex 
conjugate of rj. An account of what happens if these assumptions do not hold is given below. 
Let Wj be an eigenvector of A corresponding to the eigenvalue A,. Since Aw\= \{Wi implies 

\iWi = Awi = Awi, where the bar denotes complex conjugate, we see that if X t "=^t+i then we can 
choose Wi+i = m. It is well-known from matrix theory that because A has distinct eigenvalues any 
vector can be expressed as a linear combination of W\ , . . ., w n . So write 

v(0) = c x w\ + CiWi ■+-... + c n w„. (7) 

Taking the complex conjugate of both sides of (7) and using the fact that v{0)=v(0) we see that 
c x is real and that if Wi+i = Wi then c, + i = c,. Applying/^ ( to (7) gives 

v(t) = c x r ( wi + c 2 k! 2 W2 H- . . • + c n \l } w n . (8) 
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If A 2 is not real then X 3 must be the complex conjugate of it. That means r 2 = r 3 , 6 2 — — 3 . 
Since r\ < r 2 for i > 3,v(t) can be approximated by 

v(t) = Cir f wi + c 2 \ l 2 w 2 + c 3 \^w 3 

= Cir f Wi + c 2 r| (cos t0 2 + i sin td 2 )w 2 

(9) 
+ c 3 r|(cos £0 2 — i sin td 2 )w 3 

= C\r x W\ + r^ [cos 2 £ (c 2 u>2 + c 3 w 3 ) + sin# 2 £(c 2 w; 2 — c 3 w 3 )i] 

for sufficiently large t. 

Now C\r f Wi is, of course, the stable part of the population and the other two terms, which are 
both real by the remarks above, give rise to the fluctuations from the stable population; that is, 
they give rise to waves. It is clear from (9) that the period of the oscillation is 2ttJ \S 2 \. Also we see 
that the damping is determined by r 2 : the size of the wave relative to the total population decrease 
at a rate of r 2 \r per time period. 

If r 2 is real, then the w 3 term can be left out of the approximation in (9) and we have 8 2 = tt and 
thus the period will be two. 

The case when A is imprimitive is important for insect population. In this case r 2 = r which 
means the waves never die out. Suppose that A has exactly h eigenvalues whose absolute value is r. 
Then by Theorem 1 these eigenvalues are: 

. (h -1)6 _. {h-iW 

r, re ie , re~ ie , re i2d , re~ m , . . ., re 2 , re 2 if h is odd 

(10) 

r, re ld , re~ ,e , re l2d , re' 126 , . . ., re 2 if A is even 

where 6 = 2ir\h. The oscillations corresponding to these eigenvalues have periods of: 

A, A/2, A/3, . . .,2A/A-1 if A is odd 

(11) 
A, A/2, A/3, . . ., 2A/A = 2 if A is even. 

From this we see that the waves have periods which have harmonic ratios. Note that the entire 
sum of the waves will repeat itself every A time periods. 

We now consider what happens if we relax some of our assumptions. If there are other eigen- 
values, besides the complex conjugate of X 2 which have the same absolute value as X 2 and the 
eigenvalues of A are still distinct, then the population waves would be the sum of the waves given 
by all the eigenvalues whose absolute value is r 2 . In general the total fluctuation will be the sum 
of the waves generated by all the eigenvalues other than r. However the magnitude of oscillation 
of the eigenvalues A, with r\ < r 2 will become small as compared to the magnitude of the oscilla- 
tions of the eigenvalues with absolute value equal to r 2 . 

Now suppose r>r 2 = |X 2 | and that k 2 = \ 3 =. . . = \ k and that rj= \kj\ < \k 2 \ = r 2 for j > k. 
Let W\, . . ., w n now be a basis which gives rise to the Jordan Normal Form, see Marcus and 
Mine, (1964). That is 

Awi = rwi, Aw 2 —\ 2 w 2 ^ 

Awj = Wj-i + kjWj j= 3, . . ., k. (12) 
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As before, we expand v(0) in terms of w x , . . . , w n . Applying^ ' to v(0) and ignoring terms correspond- 
ing to eigenvalues whose absolute value is less than r 2 , we get 



v(t) = C\r t W\ + 



+ [c 2 H + c 3 t\ t 2 - 1 + 



c 4 t(t-l) w _„ 



Ck 



t(t-l). ..(t-(k-3)) ,. 



(A-2)! 



A.', 



\U)i 



+ [c 3 \' + c 4 tK- l + . . .+c k - 



Kt-l)...(t-(t-4)) ^ 
(A-3)! 
(13) 
*(*-!)• • . (t-(*-3)) 



+ . . .+ 



(A;-2)! 



r (A-2) f \ A - 2 



\|-(^-2) 



C 2 



(A: — 2)!X^" 2 



(13) 



t(t-l) . . . (t-(ifc-3)) 



■. . .-\-C k \W 2 



k 1 



*U-1). . . (f— (A — 3)) 



M>* 



Again for large t this can be approximated by 

t(t-l) . . . (t- (A -3)) 



f(£) = Cir'wi -h c/ v 



(*-2)! 



X^- (A "- 2) w;2. 



(14) 



We see from this equation that eventually the population vector does exhibit wave behavior. 
We can also see that the period of oscillation is the one which is associated with X 2 and the relative 
damping from time period t — 1 to t is 



(15) 



(t-(k-2))r 



which approaches the usual value of r 2 /r for large t. 

If there are several groups of eigenvalues all having absolute value r? then the overall wave 
will be a sum of terms as in (14), one such term for each group of eigenvalues. The behavior of these 
terms is more complicated but for large t these terms become similar to the waves obtained in (9). 

3. Stability and the Relation Between the Perron Root and Other Demographic 

Variables 

A population is said to be stable over time if the percentage of the population in each age group 
is constant over time, that is the vector of population divided by the total population is constant 
over time. 

It is well-known that a constant primitive transition matrix leads to stability over time. 

As a converse we prove the following. 
THEOREM 4: If a population exhibits stability over time then its transition matrix A(t) is constant 
if and only if the age specific birth rates are constant over time. 

PROOF: If the transition matrix A(t) is constant over time, the birth rates must, of course, be 
constant. If the age specific birth rates, aij,7=l, . . . , n, are constant, then it suffices to note 
that in a stable population the survivorship ratios aj,t_i, i = 2, . . ., n are constant. This is the case 
since stability implies each age group is increasing at the same constant rate. 

The theorem must be stated in this form since varying sets of birth rates could lead to es- 
sentially the same process if they deliver new cohorts increasing at the same rate as the population 
each time period. 
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In a one sex population the net reproduction rate (NRR) or expected number of births to a 
female allowing for the probabilities of mortality can be denned within the matrix model as the sum 

#1,1 H~ #1,2#2,1 + d\, 3#3, 2#2,1 H" • • • 

(16) 

T" fli,n#n,n-l#n-l,n-2 ' * ' a ' 2 ^' 

The rate of natural increase (RNI) of a stable population is merely the constant rate of increase 
it is experiencing in total numbers over time. Clearly RNI=\\ — 1 in a stable population, where 
• A-i = r is the maximal positive eigenvalue. (A population is said to be stationary when Xi — 1 = 0.) 

Since the NRR is somehow also a measure of how the population is increasing, intuition would 
expect a positive correlation between the NRR and RNI and thus Xi. The following theorem states 
the greatest extent to which that intuition may be formalized. 
THEOREM 5: (I) A stationary one-sex population is characterized by RNI = (Xi — 1)= (NRR— 1)= 0. 

(2) In a stable one-sex population, (i) NRR > 1 if and only ifk\ > 1 and RNI > if and only if 
NRR > X x and (ii) NRR < 1 if and only ifk,<l and RNI < if and only if NRR < k u 

PROOF: The characteristic polynomial, det (XI — A) of our matrix form (3) was computed to be 

/(X) = X n -a 1 ,iX w - 1 -ai, 2 a2,iX w - 2 -a 1 , 3 a3 5 2a2 J iX w - 3 -. . . 

. . . — ai, n &n,n-l(ln-l,n-2 • • • #21 in (4). 

It is important to note that the sum of all the subtracted coefficients is merely the NRR. 

Now, in case (1), if Xi = 1 is a root of/(X) = 0, then 1 = 1 — /(l) =NRR and conversely. 

In case (2), Ajk x has maximum eigenvalue k[ = 1 by linearity. But this means that 1 = a u i/Xj + 
«i,2#2,i/Xf + ai,3«3,2«2,i/XJ + . . . + ai, n a n , n -idn-i,n-2 • • . «2,i/X" (from the characteristic poly- 
nomial of Ajki) which is equivalent to 

Xi = #1,1 + #i,2#2,l/Xi + #i,3a3,2#2,l/Xj 4" . . .4" a l,n#n,n-ltt n _i, n _2 . . .Gte.l/X" -1 . 

(17) 
Thus Xi > 1 if and only if 

#1,1 4- #i,2#2,i 4" . . . 4" ai, n #n,n-l#n-l,n-2 • • • #2,1 > 1 

if and only if NRR > 1, and in this case NRR must be greater than k lm 

The proof is similar for Xi < 1 which completes the theorem. 

Since the same NRR could arise from many different sets of #/,/s there is no more precise 
relation between the NRR and X x . 
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